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Theoretical prediction of antilocalization of light in ultracold atomic gas samples, in the weak 
localization regime, is reported. Calculations and Monte-Carlo simulations show that, for selected 
spectral ranges in the vicinity of atomic *^Rb hyperfine transitions, quantum coherence in optical 
transitions through nondegenerate hyperfine levels in multiple light scattering generates destructive 
interference in otherwise reciprocal scattering paths. This effect leads to enhancement factors less 
than unity in a coherent backscattering geometry, and suggests the possibility of enhanced diffusion 
of light in ultracold atomic vapors. 
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I. INTRODUCTION 

An area of developing interest in fundamental atomic 
physics research is the influence of disorder on atomic 
or on mesoscopic dynamics. As an example, recent the- 
oretical results T] have indicated that disorder-induced 
transitions from a Bose-Einstein Condensate to Bose or 
to Anderson glass phases is possible, and within current 
experimental capability. Recent studies by Krug, et al. 
1^] have shown how disorder within the eigenenergy spec- 
trum of a single atom can explain anomalies in multipho- 
ton microwave ionization cross-sections. Another area 
of considerable current research interest is localization 
of light, which may be described as a disorder- induced 
phase transition in the transport properties of electro- 
magnetic radiation in strongly scattering random media. 
Investigations in this area, originally stimulated by An- 
derson ^] localization of electrons, have resulted in two 
demonstrations of light localization in condensed matter 
systems 

Recently, coherent multiple light scattering has been 
experimentally observed in samples of ultracold *^Rb 
and Sr atoms 0, Q • Corresponding theoretical results 
0, H, 0, 01 have revealed fundamental mechanisms re- 
sponsible for the observations, and numerical simulations 
have had considerable success in describing the experi- 
mental results. In all cases, the essential physical mecha- 
nisms are due to interferences in multiple wave scattering 
from the constituent atoms of the medium. For ultracold 
atoms, and under otherwise not stringent conditions, the 
interferences survive configuration averaging, thus gener- 
ating macroscopic observables. For radiation incident on 
a disordered system, the effect manifests itself as a spa- 
tially narrow ('~1 mrad) cusp-shaped intensity enhance- 
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ment in the nearly backwards direction ; this is 

referred to as the coherent backscattering (CBS) cone. 
In the present case, the disorder is principally due to the 
random locations of the atoms as found in a typical ul- 
tracold atomic gas sample. The detailed angular profile 
and magnitude of the enhancement depends on the po- 
larization of the incident and the detected light, and on 
the physical size of the sample under study. For classical 
radiation scattering from a ^ Sq —> ^ Pi atomic transition, 
the largest possible interferometric enhancement is to in- 
crease the intensity by a factor of two 0, 0, 0| over the 
incoherently scattered background light. 

Earlier we have shown that energetically remote 
atomic resonances can influence the spatial distribution 
of backscattered intensity in the cone region. For near- 
resonance scattering, these transitions appear as a spec- 
tral asymmetry, in certain polarization channels, in the 
interferometric enhancement. In addition, for CBS from 
a completely oriented atomic gas, we have found, in a 
double scattering limit, a near-resonance enhancement 
approaching the classical limit of two '6j. In the present 
paper, we present a surprising theoretical result which 
is due to the presence of off-resonance transitions. We 
find, for non resonant scattering, but over a quite wide 
spectral range, destructive interference in certain polar- 
ization channels of the coherently backscattered light. 
This results in so-called antilocalization in the op- 
tical regime, with a persistent enhancement less than 
unity. Such effects in electron transport have been widely 
studied in a range of physica l sy stems where spin-orbit 
interaction is important [H Eg. Theoretical study of 
polarization effects in coherent backscattering of s = 1/2 
massive particles by Gorodnichev, et al. ^3 have shown 
related novel effects 

In the following sections we present a brief overview 
of the physical system, including interferences generated 
by off resonance light scattering from multilevel atoms. 
This is followed by a summary of the essential theoreti- 
cal ideas, and our approach to simulate coherent multiple 
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scattering in an ultracold atomic gas. We then present 
theoretical results, in a double scattering limit, illustrat- 
ing the antilocalization effect and its physical origin. 



II. OVERVIEW OF THE PHYSICAL SYSTEM 

A. Interference in nonresonant atomic light 
scattering 

The amplitude for quasielastic scattering of light which 
is detuned several natural widths from any atomic reso- 
nance has significant contributions from all energetically 
nearby atomic transitions. The scattering cross-section 
is essentially given by the Kramers-Heisenberg equa- 
tion. To illustrate this point, consider light scattering 
from a single ground level, and two non degenerate ex- 
cited levels. As the dispersive part of the scattering am- 
plitude changes sign when tuning through a resonance, 
there will be a minimum (near-zero) in the scattering 
amp litude in the spectral region between the two levels 
[l9ll20l| . The scattering amplitude changes sign again as 
the incident radiation is tuned through this minimum. 
When we consider the complexity of real atomic transi- 
tions, the spectral location of the minimum depends on 
the polarization of the incident radiation, and on the mul- 
tipole distribution within the ground state of the scatter- 
ing atom. This behavior of the dispersion, when selected 
in one of two otherwise reciprocal multiple scattering tra- 
jectories is what is directly responsible for the existence 
of antilocalization in atomic CBS. 



B. Brief overview of the theoretical treatment 

A theory of the CBS process in an ultracold atomic gas 
has been developed recently by several groups |E BlS- 
The theoretical development essentially maintains earlier 
conceptions of weak localization in the atomic scattering 
problem [23|, and takes into account the influence of the 
optical depth and sample size on the character of the 
CBS cone. Our treatment also accounts for spec- 
tral variations on the CBS cone, and on the influence 
of far-off-resonance transitions. More recently, we have 
generalized the theoretical treatment to include effects 
associated with an angular-momentum-polarized atomic 
medium as discussed in liol . In both our earlier 
theoretical approach and in the present report, the 
general analytical development was realized by a Monte- 
Carlo simulation of coherent multiple scattering in an ul- 
tracold (T < 50 /iif ), Gaussian-shaped sample of gaseous 
^^Rb atoms confined to a magneto optical trap. These 
simulations, which included experimental parameters to 
characterize the sample, were in excellent agreement with 
our experimental results @ . 

The most important modification of the simulation 
procedure, comparing with that which we made earlier. 



is inclusion of mesoscopically averaged polarization ef- 
fects for light propagating in the atomic sample. The 
problem of light propagation through a sample consist- 
ing of atoms with arbitrary polarization in their angular 
momentum was discussed first by Cohen- Tannoudji and 
Laloe In those papers, such well known effects of op- 
tical anisotropy of an atomic vapor as birefringence, gy- 
rotropy and dichroism were connected with the formalism 
of the irreducible tensor components of atomic polariza- 
tion, i.e. with the orientation vector and the alignment 
tensor. In the Green's function formalism, the problem 
of light propagation through the polarized atomic vapor 
was later discussed in Ref. 23] . The technique devel- 
oped there allows one to obtain an analytical solution 
for the Green's function in many practically important 
cases. Basically, the retarded-type Green's function can 
be found as a solution of a certain type of Dyson equa- 
tion, which reveals the normal vacuum wave equation 
modified by the contribution of the polarization opera- 
tors on the right-hand side. These operators include all 
the polarization sensitive effects for forward propagating 
light. In turn, that lets us make a self-consistent Monte- 
Carlo simulation of the CBS process, including all the 
important polarization-sensitive effects. For more details 
we address the reader to ^ . 

Finally we emphasize that in all cases the simulations 
are made for conditions quite close to those in the ex- 
periment. These conditions include sample size, temper- 
ature, shape and density, and the characteristic intensity 
of the CBS laser beam. These conditions are such that 
recurrent scattering may be neglected, and such that sat- 
uration of the atomic transition is also negligible. In sim- 
ulations of thermal effects, and of the influence of atomic 
magnetization on the coherent backscattering enhance- 
ment, more severe conditions are used in order to illus- 
trate the possible range of influence of these effects. 



III. RESULTS 

A. Spectral variation of the CBS enhancement and 
antilocalization 

In this communication, we consider the helicity of the 
incident radiation to be antiparallel to the orientation of 
the atomic vapor. This geometry requires special prepa- 
ration since, because of optical pumping, there is a ten- 
dency to reorient the collective spin vector of the atomic 
ensemble along the beam, especially following a long in- 
teraction with the probe light, and after accumulation of 
enough Raman-type transitions. However, at this stage 
we can neglect the optical pumping mechanism and as- 
sume that most of the atoms populate the \FQ,m = —Fq) 
Zeeman state. In addition, we assume that the ensem- 
ble is located in an external magnetic field directed along 
the light beam. Thus the photons scattered via Raman 
channels are generally Zeeman frequency shifted. We as- 
sume this splitting to be quite small, and comparable or 
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FIG. 1: Diagram explaining the CBS phenomenon in the he- 
licity preserving scattering channel for double scattering of 
circular polarized light on an ensemble of *^Rb atoms ori- 
ented opposite to the helicity vector of the probe beam. In 
this example, double scattering is a combination of Rayleigh- 
and Raman-type transitions. The solid and dashed lines in- 
dicate the interfering direct and reciprocal scattering paths 
for probing between Fo — 3 ^ F — 4 and Fq = 3 ^ F = 3 
hyperfine transitions. 



less than the natural line width of the atomic transitions. 
However, because of typically very slow relaxation in the 
ground state, such splitting is enough to consider any 
Raman-type transition as an inelastic channel in the scat- 
tering process. In the following section we will explain 
how such quasi-elastic scattered modes can be resolved 
and selected in experiment. 

In Fig. n]we show double backscattering of a positive 
helicity incoming photon on a system consisting of two 
^^Rb atoms; the exit channel consists of detection of light 
also of positive helicity. The two interfering channels, 
which are shown here, repopulate atoms via Raman tran- 
sitions from the \Fq, m = —Fq) to the \Fo, m = —Fq + 2) 
Zeeman sublevel. In the direct path the scattering con- 
sists of a sequence of Rayleigh-type scattering in the first 
step and of Raman-type scattering in the second. In 
the reciprocal path, Raman- type scattering occurs first, 
and the positive helicity photon undergoes Rayleigh-type 
scattering in the second step. Since identical helicities in 
incoming and outgoing channels have opposite polariza- 
tions with respect to the laboratory frame, there is an 
important difference in transition amplitudes associated 
with Rayleigh process for these two interfering channels. 
Indeed, in the direct path the <t+ mode is coupled with 
i^o = 3^F = 4, i^o = 3^F = 3andFo = 3^i^ = 2 
hyperfine transitions. But in the reciprocal path, the cr_ 
mode can be coupled only with the Fq = 3 ^ F — 4 
transition. As we see from the diagrams shown in Fig. 1, 
where the probe light frequency ujl is scanned, for exam- 
ple, between the Fq = 3 ^ F = 4 and i^o = 3 ^ F = 3 




FIG. 2: Same as in FigQ but the double scattering is now 
combination of successive Raman-type transitions. 



transitions, a unique spectral feature is found when the 
scattering amplitudes connected the direct and recipro- 
cal scattering channels are equal in absolute value but 
have phases shifted by an angle close to tt. From an elec- 
trodynamic point of view, such conditions are realized 
when, due to the asymmetry in the Rayleigh-type tran- 
sitions, the real part of the susceptibility of the sample 
is positive for the (t_ mode and is negative for the tT+ 
mode. Since, in a first approximation, the amplitudes of 
the processes shown in Fig. 1 have opposite signs they 
will interfere destructively with anti- enhancement of the 
light scattered in the backward direction in the helicity 
preserving channel. Such an unusual behavior in the CBS 
process is connected with the Raman nature of the helic- 
ity preserving scattering channel. We reiterate that, for 
this effect to be observed, and not washed out by other 
competitive and constructively interfering channels, both 
special light polarization channels and selection of certain 
spectral domains, are required. 

In our example, as shown in Fig. 2, there are two com- 
peting channels of double Raman-type scattering which 
can interfere constructively. For equidistantly split Zee- 
man sublevels these processes contribute in the helicity 
preserving channel of the backscattering with the same 
frequency shift for the outgoing photon as for the dia- 
grams shown in Fig. ^ In an experiment, possible se- 
lection of either destructively or constructively interfer- 
ing channels can be done by focussing attention to their 
angular dependence with respect to the location of the 
atomic scatterers. For an optically thin sample the pro- 
cesses of Fig. 1 dominates for atoms located preferably 
along the probe beam direction. The respective angular 
factor is proportional to the probability to initiate the CT-i- 
or (T_ transition on the second atoms by the photon also 
emitted on the a-i- or (t_ transition by the first atom. 
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This probability is given by 

P++(^?)-P— (0)cxi(cos2 + 1)2. (1) 

In turn the processes shown in Fig. [21 are dominant if 
atoms are located preferably in the plane orthogonal to 
the probe beam. The respective angular factor is propor- 
tional to the probability to initiate the tt transition on the 
second atom by the photon emitted on the tt transition 
by the first atom: 

P^^{e) xsin'^e. (2) 

By comparing the angular distributions Q and |(5J, one 
can expect that, for a cigar-type atomic cloud, stretched 
along the probe beam direction and squeezed in or- 
thogonal directions, the destructively interfering chan- 
nels should give the dominant contribution. 

A more general way to select the destructively inter- 
fering channels can be achieved if the Zeeman sublevels 
are non-equidistantly split by applied fields. For heavy 
alkali atoms, as in the current example of *^Rb, this can 
be done by admixture of an external electric field. But it 
would be even more straightforward to accomplish this 
for lighter alkali atoms, where the splitting is always non- 
equidistant because of quadratic Zeeman effect. In this 
case the contribution of the processes shown in Fig. |21can 
be cancelled out and in double scattering the destructive 
interference can drop the output intensity down to near 
the zero level. 

However, as can be verified, the effect of destructive in- 
terference can be perfectly selected only for double scat- 
tering. The destructive interference associated with the 
diagram of Fig. 1 can be constructive for more general 
diagrams, including higher orders of scattering. We can 
expect that for a dense sample, because of accumulation 
of multiple scattering contributions of different orders, 
the effect of the destructive interference in the backscat- 
tered light still survives but can be suppressed. Realistic 
estimation of the remaining anti-enhancement factor can 
be done only numerically. 

In Fig. I^lwe show the results of numerical simulations 
of the enhancement factor in the case of pure double scat- 
tering for a spherically symmetric dense Gaussian-type 
cloud and for the Zeeman splitting in the ground state of 
O.I7. In the graphs plotted in Fig. 3 the optical depth 
has the same value of & = 1 for each detuning laser fre- 
quency wz, or corresponding detuning A = — W43. 
The solid curve shows the contribution of the processes 
depicted in Fig. The frequency of the probe laser 
ujl scans the whole upper hyperfine manifold from the 
-Fo = 3 ground hyperfine level, and there are two evident 
off-resonance spectral domains associated with antilocal- 
ization. In these spectral ranges, the enhancement factor 
is negative. It is interesting that even in such an ideal 
situation, the anti-enhancement is less than 100%. This 
is because of the polarization dependence of the refrac- 
tive index along the light ray, which is characterized by 
the properties of the Green propagation function; this 
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FIG. 3: The enhancement factor in the helicity preserving 
channel for double scattering as a function of detuning A — 
ujl — 1^4,3 for the probe laser scanning the upper hyperfine 
manifold of ^'^Rb. The solid curve shows the contribution of 
the processes depicted in Fig^and in the dashed curve the 
constructively interfering channels shown in Fig|21were added. 
Vertical lines indicate the locations of hyperfine resonances. 



essentially modifies the results of the ideal conditions re- 
quired for complete destructive interference. The chain 
curve in Fig. |31 shows how the enhancement factor would 
change if both the destructive and constructive interfer- 
ing channels shown in ^ and |2 were taken into account. 
Comparison of these dependencies suggests to us that an- 
tilocalization would be more readily observed in a system 
with non-equidistant Zeeman sublevels. 

An essential and important point is that maximal con- 
structive interference in the CBS process is normally as- 
sociated with time reversed symmetry of the scattering 
amplitude. Indeed, the reciprocal path is described by 
an amplitude which, in certain cases, can be similar to 
a time reverse amplitude of the direct path. That is the 
situation if there are only two interfering channels and 
Rayleigh type scattering. But that is not the case for 
processes appearing in the Raman channels and assisted 
by a magnetic type interaction such as hyperfine inter- 
action of nuclear and electronic spins. The destructively 
interfering amplitudes shown in Fig. 1 are not reversed 
in time. The complete time symmetry needs that the 
sign of the field to be changed for the time reverse am- 
plitude, leading to a sign change in angular momentum 
projectors. In our case the transition amplitude is in fact 
determined by the matrix elements of only the operators 
acting on electronic degrees of freedom and nuclear spin 
plays a role similar to an external magnetic field with re- 
spect to such an amplitude. In turn this means that the 
processes discussed for isolated atoms will not be seen 
without internal magnetic interactions us as a hyperfine 
interaction, which is in fact obvious as far as the ground 
state of atom would be a spin one half state in this case. 

Finally, we point out an interesting possibility sug- 
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gested by our results. It is well known 0, ^| that con- 
structive interference in recurrent scattering can reduce 
the effective diffusion constant in a randomly scattering 
medium. Alternately, this means that optical diffusive 
transport through the medium suffers from an enhanced 
resistance to the diffusive current. In direct analogy with 
our results, destructive interferences in recurrent scat- 
tering should lead to modified diffusion of light in the 
medium. Although the regime of recurrent scattering 
has not yet been reached in experiments on ultracold 
atomic gases, it is well within the current capabilities of 
experimental ultracold atomic physics. A useful observ- 
able in that case would be the analogue of the recent 
experiments of Labeyrie, et al. 24], where the time de- 
pendence of light emerging from an ultracold cloud of 
*^Rb was measured. The rate of emergence of diffusive 
light should be enhanced in the antilocalization spectral 
regime. 



B. Proposed experimental investigation using 
light-beating spectroscopy 



In the previous section it was indicated that the calcu- 
lations were done in a small static magnetic field on the 
order of a few gauss. This produces a slight ground state 
Zeeman shift which may be exploited to select different 
quasi-elastic Raman channels contributing to the multi- 
ple scattering process. The essential tool is the light beat- 
ing method, which is a common technique of high reso- 
lution spectroscopy. In an application of this method to 
observation of antilocalization, the scattered light should 
be mixed with a local oscillator in a balanced heterodyne 
detector. This coherent and high quality stable mode 
could be a portion of the monochromatic probe light. 
Then any Raman scattering channel can be observed in 
the resultant photocurrent spectrum because of the rela- 
tively low frequency of the beat note associated with the 
Zeeman shift compared with the local oscillator. With 
this approach, the Raman scattering could be observed 
as a resonant feature in the photocurrent spectrum at 
the reference frequency equal to the Zeeman splitting be- 
tween sublevels |Fo, m = —Fq + 2) and |Fo, m = —Fq). 

Here we mention only that the photocurrent spec- 
trum reproduces the spectral profile of the scattered 
light emerging the sample and responding to monochro- 
matic incident light. While probing the sample in an off- 
resonant domain, the spectral distribution is described 
by the effect of any Raman shift and also by broadening 
of the originally monochromatic coherent wave as a result 
of atomic motion. For the sake of simplicity, let us con- 
sider an optically thin atomic sample (in the frequency 
range associated with antilocalization) and ignore effects 
associated with optical depth and anisotropy; these fea- 
tures may be quite important under real experimental 
conditions. Then the contribution of single scattering is 



described by the following spectral profile 

hiu;) oc jdre'^'^"^''^'' (Jexp \2i'^{z{t) - z(0))] ^ (3) 

where wl is the probe light laser frequency and uor is 
the carrier frequency for light scattered via the Raman 
channel. The angle brackets denote averaging over the 
velocity distribution for any atom randomly moving in 
space. For steady state conditions, the stochastic func- 
tion z{t) = z(0) -|- VzT, where Vz is a random atomic ve- 
locity projected on the z-direction, and shows the atomic 
displacement along the z-dircction for a short time incre- 
ment r. The spectral profile for double scattering can be 
expressed as 



exp 



c 



^12(0)) 



(4) 

where 512(7') denotes the displacement of the scattering 
loop for any pair of atomic scatterers. This displacement 
is associated with the random motion of each atom, and 
includes all the spatial coordinates with no preference 
for any particular z-direction. The averaging in is 
extended over the velocity distribution of two atoms. In 
this expression we ignore Sagnac-type retardation effects 
for reciprocal paths; we are neglecting the small displace- 
ment of atoms during the time of light propagation be- 
tween the scatterers. 

From these equations, we see that the spectral profiles 
of single and double scattering are similar but not iden- 
tical. Thus it would be possible to resolve these terms 
in a sufficiently precise spectral analysis. The most crit- 
ical experimental aspect would be to resolve the Raman 
components in comparison with the Doppler broaden- 
ing, which should then be much less than Zeeman split- 
ting. This means that the atoms should be quite cold and 
heating effects associated with the probe light interaction 
would be undesirable. Such conditions are readily achiev- 
able using techniques of ultracold atomic physics; similar 
conditions applied to our previous measurements of co- 
herent backscattering from ultracold atomic Rb Q As 
a preliminary estimate, we compared the total intensities 
given by the integrals in ijJJ and Q and of similar expres- 
sions for higher orders of multiple scattering over uj. This 
gives us the enhancement factor for such a Raman-type 
scattering channel, which could be expressed in terms of 
cross sections, as 



X 



da- 



EF 



da^ 



da^ 



daS 



da^ 



(5) 



where the denominator contains the contributions of 
single scattering and ladder terms of all scattering or- 
ders and the numerator has an additional interference 
term. For a sample where the spectrally-dependent op- 
tical depth b in the range of the enhancement minimum 
is varied near unity, the enhancement factor changes in 
the interval of a few percent and has a rather stable min- 
imum value near —0.03. An effect of such size is readily 
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within current experimental capability. The relatively 
small value can be explained by the fact that double 
scattering cannot fully compete with the contributions 
from single scattering. As we have verified by numer- 
ical simulations, in a dense sample with b increased to 
ten, and including higher orders of scattering, the effect 
of destructive interference appearing in double scatter- 
ing also survives in the level of a few percent, but is 
partially washed out by the losses of light via other scat- 
tering channels and by effects of optical anisotropy. 

IV. SUMMARY 

A theoretical study of spectral variations in the co- 
herent backscattering enhancement factor, for a multi- 
resonance atomic gas, has been reported. Numerical sim- 
ulations of the variations, which consider the influence of 



hyperfine interferences in multiple light scattering, and 
atomic magnetization, show that antilocalization, where 
there is destructive interference in the coherent backscat- 
tering cone, can occur in the optical regime. This, in 
turn, suggests the possibility of enhanced diffusion of 
light in an ultracold atomic gas. Finally, an experimen- 
tal approach to observation of the predicted effects, using 
light-beating spectroscopy, is proposed. 
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